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Abstract 

^ , In AA collisions fusion and percolation of colour strings is studied at fixed rapidity y. Dis- 

| tiibutiou of striugs in rapidity is obtained from the observed rapidity spectra in pp collisions. 

OO ■ For y-dependence of multiplicities in Au-Au collisions good agreement is obtained with the ex- 

. isting experimental data. Predictions for LHC energies coincide with the extrapolation of the 

I data. Agreement with the data of the transverse momentum spectra requires introduction of 

' quenching into the model. 

'PhI 1 Introduction 

D I The color string model with fusion [1], [2] and percolation [3]- [6] has produced results on mul- 

' tiplicities of secondaries which are in general agreement with the existing experimental data. 

The string fusion model predicted a strong reduction of multiplicities both at RHIC and LHC 
k> I energies. At the time of the ALICE Technical proposal of 95 [7], most of model predictions, 

; (including VENUS [8], HIJING [9], SHAKER [10] and DPM [11]) for LHC energies were more 

than 4000 charged particles at central rapidity region for central {b < 3 fm) Pb-Pb collisions 
and only the prediction of the string fusion model was much lower. Since then many of the mod- 
els have lowered their predictions introducing several mechanisms, such as the triple pomeron 
coupling in DPM [12], stronger shadowing in HIJING [13] or other modifications in VENUS [14]. 
In the parton saturation picture, predictions for the central rapidity density of charged particles 
per participant for central Pb-Pb at the LHC energy range from around 15 [15] to a lower value 
around 9 [16]. Assuming that the observed geometrical scaling for the saturation momentum 
in lepton-hadron scattering is also valid for the nucleus-nucleus scattering, the value around 
9 is also obtained [17]. In percolation of strings the obtained values 7.3 [18] and 8.6 [19] are 
not far from the above ones, as expected, given the similarities between percolation of strings 
and saturation of partons [20]. In any case, these values lie above 6.4 which is obtained by 
extrapolation to the LHC energy of the values experimentally found at ^/s = 19.4, 62.8, 130 and 
200 GeV [21]. The values obtained in percolation have some uncertainty due to simplifications 
done in the calculations, mainly related to the dependence of the multiplicity on the energy of a 
simple string and on the rapidities of fused strings. One of the goals of this paper is to take into 
account both these dependencies, using as an input the rapidity and energy dependence of mul- 
tiplicities in pp collisions. We try to answer whether the charge particle density per participant 
is compatible with the one obtained in the mentioned extrapolation. We assume that strings 
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occupy different regions in the available rapidity interval. Then at different rapidities one will 
see different number of overlapping strings, depending on the rapidities of the string ends. As a 
result the string density and the process of percolation become dependent on rapidity together 
with all the following predictions for observable quantities. The central part of our derivation 
is the calculation of the number of strings at a fixed rapidity which follows from the known 
parton distributions in the projectile and target. To this aim we shall use the simple parton 
distributions employed in the review [11]. Also, we pay attention to the fragmentation rapidity 
region. As an input we take the observed limiting fragmentation scaling in pp collisions. We 
then obtain a similar scaling for A-A collisions, in agreement with the experimental data [22]. 

One of the most interesting features of the RHIC data is the suppression of high transverse 
momenta. The nuclear modification factor defined as the ratio between inclusive A-A cross 
section normalized to the number of collisions and the inclusive proton-proton cross-section is 
found to lie below unity, in disagreement with the perturbative QCD expectations. 

In the previous paper [6], in the framework of percolation of strings, a reasonable agreement 
with the data was obtained, describing A-A collisions as an exchange of clusters of overlapping 
strings. In percolation each cluster behaves like a new string with a larger tension, its value 
depending on the number of strings fused into the cluster and the cluster's transverse area. 
Fragmentation of each cluster was assumed to give rise to an exponential distribution in p^. 
Superposition of different exponential distributions then builds up a power-like distribution 
~ p^'^ with K inversely related to the magnitude of the dispersion in the number of different 
clusters. At low density, there is no overlapping of strings, thus no fluctuations and k is large. As 
the density increases, so does the string overlapping and more clusters are formed with different 
number of strings. So the dispersion increases and k decreases. Finally, at very high density, 
above the percolation threshold, there remains a single large cluster of nearly all the strings. 
Therefore, there are no fluctuations and k increases again becoming large. In this way, the 
suppression of high pT at large density follows as a result of formation of large clusters of color 
strings. In [6] we assumed that the spectrum of a simple string was a single exponential in p'^. 
However, at low string density, as in pp collisions, when fusion of strings is insignificant, the 
experimental data clearly show a power-like tail for the pT distribution. 

In this paper we study this point with more attention. Instead of the exponential distribu- 
tion we take the standard power- like parameterization for pp collisions [23]. The resulting pT 
distribution for central A-A collisions is again found suppressed at high pT but not enough to 
agree with the data. In order to describe the data we need additional suppression which would 
physically correspond to the fact that the produced particle, passing through a large cluster and 
interacting with the strong chromoelectric field, looses a part of its energy. This result is not 
unexpected. In fact a version of HIJING [24], with a string junction and doubling the string 
tension to simulate stronger color-fields, is able to explain the difference between baryons and 
mesons in the low and mid pr range but at high px some jet quenching mechanism is needed. In 
our framework quenching at high pT may be introduced in a simple phenomenological manner 
by taking the average p'^ of a cluster of n strings to grow with n more slowly than ^/n, as 
predicted for a single small cluster in absence of others. Choosing an appropriate n dependence 
of the average pj^ for clusters at a given string density allows to obtain a reasonable agreement 
with the experimental data on pT dependence in A-A collisions. 

2 pp collisions 

2.1 Multiplicities and numbers of strings 

The starting point for the calculation of fusion and percolation of strings in heavy -ion collisions 
is the distribution of strings in proton-proton collisions, where effects of fusion and percolation 
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are very small. Our strategy will be to extensively use the existing experimental data for the 
multiplicity per unit rapidity dijPP{y)/dy = ljP^{y) in pp collisions to extract the necessary 
distribution of strings in y from them. 

We recall that in the original DPM model without string fusion the multiplicity is given by a 
sum of contribution from the strings formed in the collision. In particular, in a configuration with 
n = 2k strings formed, which corresponds to the exchange of k pomerons [11], the multiplicity 
is given by 




We work in the cm system of colliding protons; Y is the overall rapidity admissible for nearly 
massless quarks. It is related to the beam rapidity as 

Y = Ybeam + ln—, (2) 

where m is the nucleon mass and /j, the quark average transverse mass. The strings are enu- 
merated according to their flavour content, that is according to which quark they are attached. 
Number 1 corresponds to the quark-diquark (qd) string, number 2 to diquark-quark (dq) string 
and all the rest correspond to sea quarks, which include ss and ss strings. Ends of strings in 
rapidity are denoted by Ui in the projectile and Wi in the target. Distributions p{ui, .■■,Un) and 
t(wi, ...Wn) give the probability to find the relevant quarks in the projectile and target proton 
respectively. Note that in the assumed notation p{u\, ...Un) gives the probability to find in the 
proton the valence quark at rapidity ui the diquark at rapidity U2 and the sea quarks at ra- 
pidities us,...Un- Distribution t(wi, ...Wn), on the other hand, gives the probability to find the 
diquark at rapidity wi, the quark at rapidity W2 and sea quarks at rapidities W3, ...Wn- Function 
Hj (y, u, w) gives the multiplicity per unit rapidity from the jth string at rapidity y provided its 
ends arc at u in the projectile and w in the target. Obviously this probability is zero if the 
string lies outside rapidity y. So iij{y,u,w) has a form 

IJ.j{y,u,w) = p{y,u,w)fij{y,u,w), (3) 

where 



p{y, u, w) = 6{u-w - yo)6{u - y)9{y - w) + 6{w -u- yo)6{w - y)9{y - u). (4) 

Here the two terms correspond to the two possibilities of the higher rapidity end of the string to 
lie on the projectile or on the target parton. The rapidity interval yo corresponds to the minimal 
extension of the string in rapidity. We take yo = 2. 

The string density dNPP(y)/dy = NPP{y) at a given rapidity is given by an expression similar 
to (1) but without jl: 

/■y/2 " " 

^n^{y)= Y[dUidWip{ui,...,Un)t{wi,...,Wn)^p{y,Uj,Wj). (5) 

•'"^/^ i=l j=l 

To analyse fusion probabilities in nuclear collisions we need to know the latter quantity. 

In principle, knowledge of the distributions p{ui, ■■■,Un) and t{wi, ...,Wn) and of string lumi- 
nosities jlj{y,u,w) allows to calculate both pf^iy) and NPP{y). This was done in the extensive 
calculations within the original DPM model [11]. However these input quantities are in fact 
poorly known and our idea is to directly relate jU^^(y) and N^^^y) using the experimental data 
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for the former. To do this we assume that string luminosities are approximately y-independent 
and the same for all type of strings: 

jlj{y,u,w) = jiQ. (6) 

This approximation has been widely used in analytical studies of string fusion. It can be justified 
for relatively long strings far from their ends, when particle production can be well described by 
the Schwinger mechanism of pair creation in a strong field. With strings of finite dimension it 
may be considered as a sort of averaging over their length and rapidity of emission. With this 
approximation we obtain a simple and direct relation between the multiplicity and number of 
strings per unit rapidity: 

lC{y) = l^,Nri^{v)- (7) 

In fact this relation for the central region was extensively used in earlier studies of string perco- 
lation. 

Relation (7) allows to find only the total number of strings per unit rapidity from the exper- 
imental data on multiplicities. However we need something more. In nucleus-nucleus collisions 
separately enter multiplicities coming from the valence strings and sea strings. Obviously one 
cannot find each of them from the experimental data. So we choose to calculate the contribution 
of valence strings from the theoretical formulas (5) and (7) and then, subtracting this contribu- 
tion from the experimental multiplicities, find the contribution from sea strings. This procedure 
can be justified by the fact that distributions of valence quarks are much better known and less 
dependent on the overall energy than the sea contribution. 



2.2 Total and sea strings from the experimental data 

We calculate the number of quark-diquark strings from (5) as 

fY/2 

Nl\y) = / dudwq^^ (u)^*) {w)p{y, u, w) 

J-Y/2 

fY/2 i-wi 

= J^ duJ^^^^dw[q^P\u)4\w)+4\u)q^P\w)), w, = mm{y,u - yo}. (8) 

Here qn^ {u) and dn^ {w) are inclusive probabilities to find a valence quark in the projectile and 
a diquark in the target at rapidities u and w respectively in a configuration with n strings. The 
second term in (8) corresponds to inverse strings whose upper ends lie on the target diquark. 
In our symmetric case the number of diquark-quark strings is obviously the same, so that the 
total number of valence strings is just twice the expression (8). 

The final number of valence strings at given y is obtained after averaging over the number 
of formed strings: 

N^y) = 2 y: ^kNgk)iy) ^ ^i^in^y))- (9) 

fc=i 

Here cu^ is the probability for the exchange of k pomerons, given by [11] 

cok = (10) 

where is the cross-section for k inelastic collisions. It is standardly taken in the K.A.Ter- 

Martirosyan model [25] 

an{s) = 2tt r bdbe-'^^^^^, (11) 
Jo nl 
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where the eikonal x corresponds to the single pomeron exchange 

x(s,6)=C(s)e-''>oW, (12) 

with ^ 

bl{s) = AR% + 4a'(lns-i^y C{s) = ^^{se-'^/^y~\ (13) 

a and a' are the pomeron intercept and slope, g is its coupling to the proton and Rn the proton 
radius. Some improvement of these ak to include the triple pomeron interaction and diffractive 
states may be found in [11]. 

To calculate the number of valence strings per unit rapidity we have to know the inclu- 
sive distributions of quarks and diquarks. Following [11] we choose the exclusive distribution 
p{ui,U2, ■■■Un) for a projectile in a factorized form 

n n 

Pn{ui,U2,...Un) = 0^(5(1 - ^Xi)Y[x1\ (14) 

i=l i=l 

where for the quark fii = 1/2 and for the diquark ^2 = 5/2 For the sea quarks and antiquarks 
we take = 1/| Inxd, where in accordance with [11] Xc = md^fs with nic = 0.1 GeV is a cutoff 
at small x. Scaling variables are related to rapidities as 

X = e-^/2+«, (15) 

Note that the distributions Pn{ui,U2, ■■■Un) are defined and normalized in the interval < x < 1, 
that is for —00 < u < Y/2. The actual strings are formed only in the part of this interval with 
u > —Y/2. This circumstance is inessential for valence quarks whose distributions rapidly vanish 
towards small values of x. For the distributions in the target one has only to invert the rapidities 
u — ^ —w in (15). Integration over the scaling variables of unobserved partons gives the desired 
inclusive distributions. For the valence quark we find 
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For the diquark 



f W-c.^ U ^) > 0Fr(5/2 + (n-2)/x)- ^^^^ 



d{n)( ^ ^ cwx5/2n _ a.^-l/2+(n-2)/. ^ 4r(3 + (n - 2)fi 

dp [X) CdX X) , Cd 3^r(l/2 + (n - 2);Lx) ■ ^ 

After the averaged valence string number is found according to (9) we have to transform it 
into the valence multiplicity using (7). The value of /xq can be found from the observed plateau 
hight assuming that at y = all strings contribute. Their average number can be found from (9) 
as = 2{k). As a result we find values of ^0 slowly rising with energy and visibly saturating 
at TeV energies. In Fig. 1 we show these values extracted from the data [23] together with 
their extrapolation to the LHC energies in the assumption that the plateau in the pp multiplicity 
distribution rises linearly with Ins. The obtained multiplicities from valence strings vanish in 
the fragmentation region too slowly as compared to tlie experimental data, and at y — '^beam 

> 

become greater than the latter. This is obviously related to our assumption of a constant string 
luminosity throughout the string length, whereas it should go to zero at its ends. Put in other 
words, in our approach the total energy is conserved in its division between different strings 
due to the (5- function in (14) but it is not conserved inside each separate string, since near the 
string end its luminosity should vanish. To cure this defect in a simple manner we just assume 
that as soon as the calculated valence contribution becomes larger that the data we substitute 
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Figure 1: String luminosity as a function of energy 



the former by the latter, assuming that in this deep fragmentation region sea strings do not 
contribute at all. 

With thus obtained valence contribution to the multiplicities we find the sea contribution 
just as the difference between the total and valence one. Dividing it by jUo we find the number 

of sea strings per unit rapidity in pp collisions N^''f'f{y). In fact wc need not exactly this number 
but the one in the assumption that all strings are of the sea type, which is obtained from it by 
rescaling 
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NPP - 2 

This quantity is shown in Fig. 2 for different values of the collision energy -s/s. 



(18) 



3 hA and AA collisions 

Generalization to hA and AA collisions is straightforward and follows [11]. At fixed impact 
parameter b one introduces the average numbers of participants 2i/^(6) and collisions i'{b). Then 
the number of strings in AA collisions at given b and y is 

NAAib, y) = Vpar{b)N^{y) + {vcolih) - I^par{b))N%y). (19) 

Here N^^i^y) is obtained from the observed multiplicity in pp collisions according to (7) and 
N^{y) is given by (18). 

One can easily further generalize (19) to collisions of different nuclei (see [11]). 

In the Glauber approach the numbers Vpar and Vcoi for AA collisions are obtained as follows 

/ (fb'TAih') ( 1 - e-^'^'^^^^-^'A a (b) 
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Figure 2: Sea string contribution to multiplicities per nucleon as a function of rapidity. Curves 
from bottom to top correspond to energies 19.4, 62.8, 130, 200 and 6000 GeV. 



where a is the total pp-cross-section, Ta(6) is the nuclear profile function normalized to unity 
and 

TAA{h) = J d%'TA{h')TA{h - h'). (21) 
For hA collision, as mentioned, Upar = 1 and 



aTA{b) 



1 _ (,-AaTA{b) ■ 



(22) 



Note that for AA collisions in the above formulas the denominator is written in the so-called 
optical approximation [26]. As is well-known, it works reasonably well except close to the 
nucleus boundary, where the collision numbers obtained from (20) may be quite deceptive. 

It is customary to take the profile function r^(6) generated by the Woods-Saxon nuclear 
density. However for our purpose it is more convenient to assume the nucleus to have a well- 
defined radius, which allows to determine the interaction area in the transverse plane as just the 
area of the overlap. For this reason we take the nucleus as a sphere of radius Ra = A^/^ ■ 1.2 
fm, which gives 



TA{h) = e{RA-h)- 



Va 



(23) 



where Va is the nuclear volume. With this choice the most peripheral collisions occur at 6 = 2Ra 

with Vcol = Vpar = 1- 



4 Multiplicities and pr distributions 

In our previous studies of string fusion we always stressed that it can only occur in the common 
rapidity interval. However our attention was mostly centered on the central rapidity region 
where all (or nearly all) strings contribute, so that the requirement of common rapidity interval 



8 



was of no relevance and strings could be considered as of practically infinite length in rapidity. 
Now we study the fusion process in more detail. At a fixed rapidity y only strings which pass 
through this rapidity can fuse. Formulas of the previous sections allow to find the original 
number of strings N{h, y) stretched between the projectile and target at fixed rapidity layer y 
and impact parameter b. According to the percolation colour strings scenario these strings in fact 
fuse into strings with higher colour. The intensity of fusion is determined by the dimensionless 
percolation parameter 77 proportional to the string density in the interaction area 

where sq = t^i^o is the transverse area of the string and S{b) is the interaction area, that is, the 
overlap area in case of AB collisions. Obviously in our case the percolation parameter depends 
both on b and y. At r/ ~ 1.2 -j- 1.3 fusion of strings leads to their percolation and formation 
of macroscopic string clusters. This phenomenon will take part only in restricted intervals of b 
and y, predominantly at central collisions and rapidities, where the effects of string fusion and 
percolation will be most noticeable 

Considering the case of AA collisions at reasonably high energies we shall assume the total 
number of strings high enough to allow use of the thermodynamic limit, in which the total areas 
of n-fold fused strings Sn become distributed according to the Poisson law with (n) = r){b,y): 

5„(6,y) = 5(6)e-''('''-)^^^^. (25) 

Due to averaging of the direction of colour, the n-fold fused string emits the number of particles 
which is only ^/n times greater that the simple string. So the total production rate at fixed b 
and y will be given by 

MM) = Mo^e-(^-)i:V^^, (26) 

where jUo is the production rate from the single string, which, as stated above, we assume to be 
independent of y but dependent on energy. 

As to the pT distribution, we use a slightly generalized model introduced and discussed in 
[27], in which the normalized probability Wn{p) to find a particle with transverse momentum p 
emitted from the n-fold fused string is given by 

^nip)= ^^--l^^7~^\ -^T- (27) 

2-iTpi ^p + pj 

Here for n = 1 the parameters are determined by the experimental data on pp collisions: 

Pi = 2 GeV/c, Ki = 19.7 - 0.86 In E^rn (28) 

and Ecm is the cm. energy in GeV. With n > 1 from the string fusion scenario it follows that 
the average transverse momentum squared of the particles emitted from the n-fold fused string 
is "n}^'^ greater than for a single string. This gives a relation between p^ and Kn 

p2 _ 1/2 2 (fcn-3)(fcn-4) 

So the distribution from n-fold string is fully determined by the n-dependence of In [27] 
the simplest choice of n-independent «;„ was used. However this simple choice does not allow to 
obtain the pT dependence in agreement with the data at RHIC. To improve our description, we 
introduce corrections to the original string picture which correspond to non-linear phenomena 
in string clustering and influence both and the behaviour of Pn 
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In fact the value of k controls the difference of the distribution from the purely exponential 
one, passing into the latter at very large k. For a single string a finite ki may be thought 
of as a result of fluctuations in the string tension (or equivalently its transverse area) [28]. 
One may expect this fluctuations to grow as many string fuse and so the value of k„ should 
fall with n. However there is another effect acting in the opposite direction. As many strings 
fuse into clusters, multiple interactions of emitted particle inside the clusters should lead to 
thermalization of the particle spectra making it closer to an exponential. Thus eventually at 
large n parameter k„ should grow to large values. Naturally we cannot determine the exact 
form of the n dependence of Kn from purely theoretical reasoning. We can only think that the 
change from fall to growth should occur in the vicinity of the percolation threshold and that in 
any case k cannot be smaller than 4 to have a convergent < >. In practice we take k„ as 

K„ = + o(n - 1) + 6(n - 1)^ (30) 

and try to adjust a and b to get a better agreement with the experimental data. In fact the 
results are not very sensitive to the choice of a and b provided they are taken to have a behaviour 
of Kn in agreement with the above general theoretical observations. However the generalization 
(30) is not sufficient to bring our predictions in agreement with the observed quenching of the 
ratios Raa at RHIC. To this aim wc have to introduce some quenching also in the string picture 
at large values of ry. It corresponds to the fact that passing through a large cluster volume 
and interacting with the strong chromoelectric field the produced particles loose a part of their 
energy [29] . On our phenomenological level it would correspond to the behaviour of the average 
transverse momentum squared as 

< / >„= n"" < p2 (31) 
with the exponent a„ less than 1/2 and diminishing with n. Similarly to (30) we parameterize 

an = ^ + c{n-l)+d{n-lf. (32) 

The comparison with the experimental data determines the optimal fit for the parameters a, b, 
c and d. The resulting values for Kn and are shown in Figs. 3 and 4. 

Averaging with the distribution in n we find the final distribution in p from the fusing strings 
at given b and y as 

wip,b,y)=a{b,y)J2Mp)^^^^^, a-Hb,y) = e^^'''y^-l (33) 
(the change in the normalization is due to the restriction n > 1). 

5 Numerical results 

We studied Au-Au collisions at energies 19,4, 62.8, 130, 200 and 6000 GeV corresponding to 
the existing experimental data and expected at LHC. Using our results on the string numbers 
in pp collisions wc calculated their numbers in nucleus-nucleus collisions at a given y. Knowing 
these numbers and also numbers of participants and collisions wc then determined values of the 
percolation parameter r]{y, b) at different rapidities and impact parameters. We have taken the 
transverse radius of the single string 0.3 fm. In Fig. 5 and 6 we illustrate values of r]{y, b) as a 
function of y for central collisions and as a function of b at mid-rapidity. As one can observe, 
at RHIC and LHC energies these values arc quite large, far beyond the percolation threshold. 
As a result one finds a very substantial reduction in multiplicity calculated according to Eq. 
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Figure 3: Parameter k„ in the distribution (27) for the n-fold fused string. Curves from top to 
bottom correspond to energies 19.4, 62.8, 130, 200 and 6000 GeV. The curves end at maximal 
n reached at these energies 

(26) with luminosities determined from pp colhsions. In Fig. 7 and 8 we show multiplicities as 
a function of y for central collisions and function of b at midrapidity. They agree rather well 
with the existing data both in form and absolute values. In Fig. 8 a rather sharp change is seen 
in the periphery of the nuclei, between b/2RA = 0.8 and unity. The RHIC data do not exhibit 
such a sharp saturation. This behaviour is a direct consequence of using in our calculations 
the nuclear profile function (23) corresponding to the step function for the nuclear density. A 
more realistic profile function would lead to a smoother transition from the periphery to center, 
although requiring a more complicated definition of the interaction area. 

Our prediction for the plateau at LHC energy (divided by fpor) is around 9 for the sum 
of charged and neutral particles (that is 6 for charged), which is lower than derived in other 
publications [13]-[17]. The plateau height for charged particles in central Au-Au collisions as 
a function of energy is illustrated in Fig. 9. Of course our results are directly related to the 
chosen string radius and go upward if it is lowered. However then we loose the agreement with 
the existing data. 

Note that our calculations also reproduce quite well the behaviour in the fragmentation 
region (limiting fragmentation). For this the dependence of the string luminosity on energy 
proved to be quite important. Without it the nice linear dependence of multiplicities in the 
fragmentation region is spoiled and the line is widened into a band. 

To clearly see the effect of string fusion in Fig 10 we show the multiplicities at 6 = without 
fusion. Their values are several times greater than with fusion and do not agree with the 
experimental data at all. 

Passing to the pr distributions in Fig. 11 we show the ratios Raa{pt) for central collisions 
at midrapidity. The curve for 200 GcV served to determine our parameters a, b, c and d in (30) 
and (32). Our predictions for the LHC energy show a behaviour similar to the RHIC energy 
with a still more pronounced quenching effect. In the fragmentation region we prefer to show 
the ratios with normalization respective to the number of participants, since in this region 
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Figure 4: Power a„ for the n dependence of the average transverse momentum (31) of the n-fold 
fused string 

the multiphcities are roughly proportional to i/par due to low string densities. Fig. 12 shows 
that these ratios are close to unity and may only fall a little below unity at LHC energies. 

6 Conclusions 

In the framework of percolation of strings we have obtained a strong reduction of multiplicities 
at LHC, much larger than the rest of models but in agreement with the extrapolation from 
the SPS and RHIC experimental data. Due to similarities between percolation of strings and 
saturation of partons, it would be interesting to explore the possibility for further reduction 
of multiplicities in the saturation approach. In order to describe the energy dependence of 
multiplicities in AA collisions we need a rather large transverse size of the elementary string 0.3 
fm. This enhances the interaction of strings and so cluster formation, which leads to stronger 
reduction of multiplicities. In our calculations we have used the standard optical approximation 
to compute the numbers of participants and collisions. This approximation enhances the number 
of collisions for peripheral collisions in comparison with Monte-Carlo evaluations and leads to 
some uncertainties also for central collisions. For this reason, our results must be regarded to 
have an uncertainty in the range of 10%-15% . So the string transverse size can be lower if the 
number of strings is in fact lower. 

Taking limiting fragmentation scaling for pp collisions as an input, we have found the same 
behaviour for AA collisions, which is confirmed experimentally up to the RHIC energies. Our 
calculations predict that limiting fragmetation scaling also remains approximately valid at the 
LHC energy (with a 5% suppression compared with to SPS or RHIC, see Fig. 7). 

We have been able to describe reasonably well the high transverse momentum spectrum 
at different energies ranging from SPS to RHIC. A large suppression is predicted for LHC. In 
order to obtain such an agreement, in addition to the usual effects of string clustering, such 
as reduction of the effective number of independent color sources and suppression of transverse 
momentum fluctuations, we need a shift of the pT spectrum due to energy loss. Considering 
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Figure 5: Percolation parameter rj in central Au-Au collisions as a function of rapidity. Curves 
from bottom to top correspond to energies 19.4, 62.8, 130, 200 and 6000 GeV. 

string fusion as an initial state effect (before particle production), a final state effect is needed 
to account for the observed suppression, similarly to jet quenching in the QCD picture. 
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Figure 7: Multiplicities (charged plus neutrals) per 1/2 number of participants in central Au-Au 
collisions as a function of rapidity. Curves from bottom to top correspond to energies 19.4, 62.8, 
130, 200 and 6000 GeV. 
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Figure 8: Multiplicities (charged plus neutrals) per 1/2 number of participants in Au-Au colli- 
sions at midrapidity as a function of impact parameter. Curves from bottom to top correspond 
to energies 19.4, 62.8, 130, 200 and 6000 GeV. 
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Figure 9: The plateau height for charged particles per 1/2 number of participants in Au-Au 
collisions at midrapidity clS Si function of energy 
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Figure 10: Multiplicities (charged plus neutrals) per 1/2 number of participants in central Au- 
Au collisions at midrapidity without fusion. Curves from bottom to top correspond to energies 
19.4, 62.8, 130, 200 and 6000 GeV. 
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Figure 11: Nuclear factor Raa for central Au-Au collisions at midrapidity. Curves from top to 
bottom correspond to energies 19.4, 62.8, 130, 200 and 6000 GeV. 
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Figure 12: Nuclear factor for central Au-Au collisions in the fragmentation region y = O.IY. 
Curves from top to bottom correspond to energies 19.4, 62.8, 130, 200 and 6000 GeV. 



